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1  Executive  Summary 

The  lack  of  a  systematic  methodology  for  the  design  of  feedback  laws  capable  of  controlling 
complex  dynamical  systems  has  been  a  limiting  factor  in  several  current  and  emerging  DoD 
missions.  The  research  carried  out  by  the  principal  investigator  in  this  three  year  research 
effort  has  focused  on  analyzing  and  computing  the  steady-state  behavior  of  controlled  com¬ 
plex  dynamical  systems.  The  control  of  such  systems  typically  involves  both  the  design  of 
feedback  laws  which  use  the  current  value  of  the  state  and  the  design  of  dynamical  systems 
which  produce  estimates  of  the  current  state. 

One  of  the  fundamental  discoveries  made  during  this  research  effort  concerns  the  steady- 
state  behavior  of  a  dynamical  system  which  provides  estimates  of  the  current  state,  propa¬ 
gating  smaller  amounts  of  data  than  required  by  the  standard  Kalman  filter.  The  analysis 
of  this  estimation  scheme  led  to  an  unanticipated  discovery  about  the  geometry  of  certain 
classes  of  linear  systems.  The  geometric  properties  which  were  discovered  in  this  way  led  to 
the  solution  of  an  outstanding  problem  in  linear  systems  theory  with  applications  in  speech 
synthesis,  voice  recognition  and  signal  processing.  These  advances  were  supported  by  com¬ 
putational  methods  developed  in  this  research  effort  and  which  are  documented  in  this  final 
report  and  in  a  patent  application. 

Another  fundamental  discovery  made  during  this  research  effort  was  focused  on  the  first 
part  of  controller  design  discussed  above,  for  a  class  of  distributed  parameter  systems.  The 
principle  results  obtained  quantify  properties  of  attractors  and  the  steady-state  behavior  of 
solutions  for  a  controlled  Burgers’  equation,  which  is  often  used  as  a  simple  model  for  tur¬ 
bulence.  The  principal  investigating  team  also  discovered  results  concerning  the  asymptotic 
behavior  of  linear  distributed  parameter  systems  undergoing  harmonic,  or  periodic,  forcing. 
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1.1 


The  Rational  Covariance  Extension  Problem  and  Speech  Syn¬ 
thesis 


In  the  course  of  our  studying  the  design  of  nonlinear  observers  we  discovered  a  description 
of  the  complete  phase  portrait  of  the  Kalman  filter,  viewed  as  a  nonlinear  dynamical  system 
on  the  space  of  positive  real  functions.  This  discovery  in  turn  lead  to  the  solution  of  an 
important  problem  in  signal  processing.  In  particular,  in  [7]  we  were  able  to  give  a  complete 
parameterization  of  all  shaping  filters  which  match  a  given  finite  window  of  correlation  coeffi¬ 
cients.  This  was  an  important  open  problem  in  signal  processing,  first  formulated  for  speech 
synthesis  (and  intensively  studied)  by  researchers  at  the  Phillips  Laboratories  in  Europe.  In 
[14,  15],  T.  T.  Georgiou  proved  that  to  each  choice  of  partial  covariance  sequence  and  nu¬ 
merator  polynomial  of  a  modeling  filter  there  exists  a  rational  covariance  extension  yielding 
a  pole  polynomial  for  the  modeling  filter,  and  he  conjectured  that  this  extension  was  unique 
so  that  it  provides  a  complete  parameterization  of  all  rational  covariance  extensions.  We 
note  that  this  problem  has  a  long  history,  with  antecedents  going  back  to  potential  theory  in 
the  work  of  Caratheodory,  Toeplitz  and  Schur  [10,  11,  32,  31],  and  continuing  in  the  work  of 
Kalman,  Georgiou,  Kimura  and  others  [19, 15,  22].  This  problem  has  been  of  more  recent  in¬ 
terest  due  to  its  significant  interface  with  problems  in  signal  processing  and  speech  processing 
[12,  9,  26,  21]  and  in  stochastic  realization  theory  and  system  identification  [2,  33,  23].  Our 
recent  work  [7],  extends  the  result  of  Georgiou  and  answers  his  conjecture  in  the  affirmative. 
This  work  has  shed  new  light  on  the  stochastic  (partial)  realization  problem  [6]  through 
the  development  of  an  associated  Riccati-type  equation,  whose  unique  positive  semi-definite 
solution  has  as  its  rank  the  minimum  dimension  of  a  stochastic  linear  realization  of  the  given 
rational  covariance  extension.  However,  our  proof  was  not  nonconstructive.In  [8],  we  have 
been  able  to  give  a  constructive  proof  of  Georgiou’s  conjecture,  which  provides  an  algorithm 
for  solving  the  problem  of  determining  the  unique  pole  polynomial  corresponding  to  the  given 
partial  covariance  sequence  and  the  desired  zeros.  In  this  work,  which  was  motivated  by  the 
effectiveness  of  interior  point  methods  for  solving  nonlinear  convex  optimization  problems, 
we  first  recast  the  fundamental  problem  as  an  optimization  problem. 

Historically,  the  only  solution  to  rational  covariance  extension  problem  for  which  there 
has  been  simple  computational  procedures  is  the  so  called  maximum  entropy  solution,  which 
is  the  particular  solution  that  maximizes  the  entropy  gain.  During  this  research  period,  we 
were  able  to  demonstrate  that  the  infinite-dimensional  optimization  problem  for  determining 
this  solution  has  a  simple  finite-dimensional  dual.  This  motivated  the  introduction  of  a 
nonlinear,  strictly  convex  functional  defined  on  a  closed  convex  set  naturally  related  to 
the  covariance  extension  problem.  We  were  able  to  show  that  any  solution  of  the  rational 
covariance  extension  problem  lies  in  the  interior  of  this  convex  set  and  that,  conversely, 
an  interior  minimum  of  this  convex  functional  will  correspond  to  the  unique  solution  of  the 
covariance  extension  problem.  Our  interest  in  this  convex  optimization  problem  is,  therefore, 
twofold:  as  a  starting  point  for  the  computation  of  an  explicit  solution,  and  as  a  means  of 
providing  an  alternative  proof  of  the  rational  covariance  extension  theorem. 

Concerning  the  existence  of  a  minimum,  we  showed  that  this  functional  is  proper  and 
bounded  below,  i.e.,  that  the  sublevel  sets  of  this  functional  are  compact.  Prom  this,  it 
follows  that  there  exists  a  minimum.  Since  uniqueness  follows  from  strict  convexity  of  the 
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functional,  the  central  issue  which  needed  to  be  addressed  in  order  to  solve  the  rational 
covariance  extension  problem  was  whether,  in  fact,  this  minimum  is  an  interior  point.  Indeed, 
our  formulation  of  the  convex  functional,  which  contains  a  barrier-like  term,  was  inspired  by 
interior  point  methods.  However,  in  contrast  to  interior  point  methods,  the  barrier  function 
we  introduced  did  not  become  infinite  on  the  boundary  of  our  closed  convex  set.  Nonetheless, 
we  were  able  to  show  that  the  gradient,  rather  than  the  value,  of  the  convex  functional  became 
infinite  on  the  boundary.  The  existence  of  an  interior  point  which  minimizes  the  functional 
then  follows  from  this  observation. 


Some  Numerical  Examples 

Given  an  arbitrary  partial  covariance  sequence  Co,  C\, . . . ,  Cn  and  an  arbitrary  zero  polyno¬ 
mial  cr(z),  the  constructive  proof  of  Georgiou’s  conjecture  provides  algorithmic  procedures 
for  computing  the  corresponding  unique  modeling  filter,  which  are  based  on  the  convex 
optimization  problem  to  minimize  the  functional  tp  :  Rn+1  -»  R,  defined  by 


<p(qo,  01>  •  •  •  .  0n)  =  Co?0  +  Cl 01  + - b  CaQn 


over  all  00)  9ii  •  •  •  5  9n  such  that 


Q(etd)  =0o  +  0i  cos  0  +  02  cos  2 6  + - b  qn  cos  nO  >  0  for  all  9. 


(1.1.1) 


(1.1.2) 


over  all  0o,  0i, . . . ,  qn  such  that  (1.1.2)  holds. 

In  general  such  procedures  will  be  based  on  the  gradient  of  the  cost  functional  ip,  which, 
is  given  by 


Qq  (00)  01)  •  •  •  )  0n)  —  Cfc  Cfc 

where 

'Ck  =  hS-f°r^Ai6  for*  =  <u'2’--n 

are  the  covariances  corresponding  to  a  process  with  spectral  density 

Q(e'<>)  =g°  +  2Z+cos(^). 


(1.1.3) 


(1.1.4) 


(1.1.5) 


The  gradient  is  thus  the  difference  between  the  given  partial  covariance  sequence  Co,  C\, . . . ,  Cn 
and  the  partial  covariance  sequence  corresponding  to  the  choice  of  variables  qo,Qi,  •  •  ■  ,qn  at 
which  the  gradient  is  calculated.  The  minimum  is  attained  when  this  difference  is  zero. 

The  following  simulations  have  been  done  by  Per  Enqvist,  using  Newton’s  method  (see, 
e.g.,  [24,  27]),  which  of  course  also  requires  computing  the  Hessian  (second-derivative  matrix) 
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in  each  iteration.  An  straight-forward  calculation  shows  that  the  Hessian  is  the  sum  of  a 
Toeplitz  and  a  Hankel  matrix.  More  precisely, 


(QOi  Qh  •  •  •  >  Qn)  —  ^i—j)  j  1,  2,  .  .  .  ,  TZ,  (1.1.6) 

where 

1  r *  |  [2 

*=-y  for  fc  =  0,1,2,  (1.1.7) 

and  d-k  =  dk ■  Moreover,  d0,  di,  d2, . . . ,  d2n  are  the  2n  +  1  first  Fourier  coefficients  of  the 
spectral  representation 

\  f  i6\\2  00 

Q(e»)*  =  d*  +  2  dk  cos (**)•  (1-L8) 

The  gradient  and  the  Hessian  can  be  determined  from  (1.1.3)  and  (1.1.6)  respectively  by 
applying  the  inverse  Levinson  algorithm  (see,  e.g.,  [28])  to  the  the  appropriate  polynomial 
spectral  factors  of  Q(z )  and  Q(z)2  respectively  and  then  solving  the  resulting  linear  equations 
for  Co,  Ci, . . . ,  Cn  and  d0,  di,  d2,..-,  d2n\  see  [13]  for  details. 

To  illustrate  the  procedure,  let  us  again  consider  the  sixth  order  spectral  envelopes  of 
Figure  1  and  2  together  with  the  corresponding  zeros  and  poles. 


Figure  1:  Spectral  envelope  of  a  maximum  entropy  solution. 
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frequency 


Figure  2:  Spectral  envelope  obtained  with  appropriate  choice  of  zeros. 

Hence  Figure  3  illustrates  the  periodogram  for  a  section  of  speech  data  together  with 
the  corresponding  sixth  order  maximum  entropy  spectrum,  which,  since  it  lacks  finite  zeros, 
becomes  rather  ’’flat”.  The  location  of  the  corresponding  poles  (marked  by  x)  in  the  unit 
circle  is  shown  next  to  it.  The  zeros  (marked  by  o)  of  course  all  lie  at  the  origin. 


frequency 


Figure  3 


Now  selecting  the  zeros  appropriately  as  indicated  to  the  right  in  Figure  4,  we  obtain 
the  poles  as  marked,  and  the  corresponding  sixth  order  modeling  filter  produces  the  spectral 
envelope  to  the  right  in  Figure  4.  We  see  that  the  second  solution  has  a  spectral  density  that 
is  less  flat  and  provides  a  better  approximation,  reflecting  the  fact  that  the  filter  is  designed 
to  have  transmission  zeros  near  the  minima  of  the  periodogram. 


Figure  4 


1.2  Structure  of  Attractor  for  a  Boundary  Controlled  Burgers’ 
Equation 

Motivated  in  part  by  problems  of  flow  control  and  combustion  control,  where  nonlinear  effects 
actually  can  improve  mixing,  there  has  been  considerable  attention  given  in  the  literature 
to  the  study  of  asymptotic  or  steady  state  properties  of  solutions  of  nonlinear  distributed 
parameter  systems,  such  as  Navier-Stokes  equations,  which  contain  both  nonlinear  convective 
terms  and  diffusive  terms.  To  this  end,  part  of  our  research  effort  concerns  the  effect  of 
boundary  control  on  the  structure  of  attractors  for  Burgers’  equation.  In  the,  soon  to  be 
published,  work  [123],  we  give  a  rigorous  mathematical  justification  of  the  fact  that  there  can 
be  multiple  stationary  points  contained  in  the  global  attractor.  Partial  numerical  concerning 
this  possibility  were  first  results  announced  in  [95].  The  analysis  in  this  paper  is  based  on  an 
extension  of  the  classical  analysis  of  Sturm-Louiville  boundary  value  problems  using  the  so- 
called  Priifer  transformations  or  polar  coordinates.  We  should  comment  that  various  partial 
results  were  obtained  earlier  using  a  variety  of  approaches  but  none  of  these  approaches 
provided  the  more  complete  picture  presented  here. 

In  our  work  [123]  we  show  that  for  a  special  class  of  forcing  terms  as  the  gain  parameters 
in  the  boundary  feedback  control  are  increased  from  small  positive  values  to  large  positive 
values,  the  number  of  stationary  solutions  vary  from  three  to  one.  This  is  in  keeping  with 
our  results  on  conergence  of  attractors,  discussed  in  the  last  subsection,  where  it  is  shown 
that  the  zero  dynamics  systems  has  a  single,  global,  asymptotically  stable,  equilibrium. 

Consider  the  controlled  viscous  Burgers’  system 

wt  -  ewxx  +  wwx  =  f(x), 
w  =  w(x,  t ),  x  E  (0, 1),  t  >  0, 

-wx(0,  t)  =  u0{t),  wx(l ,t)  =  ui(t),  (1.2.1) 

w(x,  0)  =  (j){x ), 

2/o(*)=w(M),  yi(t)  =  w(l,t), 
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where  u0(t),  U\(t)  are  boundary  inputs,  yo(t),  yi(t)  are  boundary  outputs,  f  €  L2 (0,1)  is  an 
external  forcing  term  modeling  an  unknown  disturbance. 

Formally  introducing  proportional  error  feedback,  with  ko  =  k\  =  k 

u0  =  -ky0,  ui  =  -kyu  (1.2.2) 

with  feedback  gain  k  >  0  we  obtain  the  closed  loop  Burgers’  system 

wt  -  ewxx  +  wwx  =  /, 
x  e  (0,1),  t  >  0, 

—wx(0  ,t)  +  kw(Q,t)  =  0, 
wx(l,t)  +  kw(l,t)  =  0, 
w(x,  0)  =  (f)(x ). 

In  our  work  we  are  interested  in  forcing  terms  f(x)  possessing  a  certain  symmetry 
property  and  an  additional  definiteness  property.  These  assumptions  considerably  simplify 
the  analysis  of  the  resulting  stationary  problem. 

Assumption  1.2.1.  We  assume  that  f  is  an  odd  function  about  x  =  1/2  in  the  interval 
[0, 1],  i.e., 

f(x)  =  -/(l  -  x)  for  x  e  [0, 1], 

and  we  introduce  the  terminology  ‘antisymmetric  about  1/2”  or  simply  “antisymmetric”  to 
describe  such  a  function. 

We  will  also  assume  that 


f(x)  >  0  for  x  €  [0, 1/2). 

An  important  fallout  of  the  antisymmetry  condition  is  that  it  is  preserved  by  solutions 
of  (1.2.3). 

The  stationary  Burgers’  system  associated  with  (1.2.3)  is 

-ewxx(x)  +  w(x)wx(x)  =  f(x),  (1.2.3) 

^i(0)  —  kw(  0)  =  0, 
wx(l)  +  kw(l)  =  0. 

Our  approach  proceeds  as  follows.  We  first  integrate  the  stationary  Burgers’  equation 
and  then  introduce  the  Riccati  (or  Hopf-Cole)  transformation  to  reduce  the  stationary  Burg¬ 
ers’  equation  to  a  second  order  linear  equation  containing  a  spectral  type  parameter.  Un¬ 
fortunately,  as  mentioned  above,  the  boundary  conditions  are  transformed  into  nonlinear 
boundary  conditions.  In  the  stationary  case  these  boundary  conditions  factor  into  products 
of  boundary  conditions  which  appear  to  be  of  Sturmian  type. 

Integrating  the  differential  equation  (1.2.3)  over  the  interval  [0,  x],  we  arrive  at 

-ew'(x)  +  =  f  f(s)  ds  +  c,  (1-2.4) 

2  Jo 
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(1.2.5) 


where,  c  is  a  constant  of  integration.  Using  the  boundary  conditions,  we  see  that 

c  =  —ekw(0)  +  ui(0)2/2 

is  actually  not  arbitrary.  The  parameter  c  plays  an  important  part  in  our  analysis. 

Equation  (1.2.4)  is  a  Riccati  ordinary  differential  equation  and  one  classical  approach 


to  solving  this  equation  is 
solution  in  the  form 

to  introduce  the  so  called  Riccati  transformation,  i.e., 

we  seek  a 

.  .  v'(x) 

w(x)  =  -2e— 7— . 

v(x) 

(1.2.6) 

Using  (1.2.6),  the  differential  equation  (1.2.4)  is  transformed  into 

v"(x)  -  ( F(x )  +  A)  v(x)  =  0 

(1.2.7) 

where 

(1.2.8) 

and 

A-  — 

A  2e2' 

(1.2.9) 

The  boundary  conditions  are  transformed  into  nonlinear  boundary  conditions,  which, 
after  using  the  equations  (1.2.7)  and  (1.2.8),  can  be  written  as 


-cv(0)2  +  2v'{0)2e2  +  2kv'{0)v{0)e2  =  0,  (1.2.10) 

-cv( l)2  4-  2v'(l)2e2  -  2kv' (l)v(l)e2  =  0. 

These  conditions  can  be  factored  into  products  of  conditions  that  appear  to  be  of  “Sturm- 
Liouville”  type: 

(t/(0)  -  l+v{0))  (v'{0)  -  l~v{0))  =  0  (1.2.11) 

(u'(l)  +  Z+v(l))  (w'(l)  +  Uu(l))  =  0. 

(1.2.12) 


Unfortunately,  the  parameters  depend  on  the  “spectral”  parameter  A: 

(±  =  !*(*,  A)  =  ±  ^(j)2  +  A.  (1.2.13) 

Considering  the  various  possibilities  we  arrive  at  four  systems  that  can  deliver  stationary 
solutions 


v"(x )  -  ( F(x )  +  A)  v(x)  =  0 


(1.2.14) 
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(1.2.15) 


t/(0)  -  l+v( 0)  =  0,  t/(l)  +  /+u(l)  =  0, 

v'(0)  -  l+v(0)  =  0,  v'(l)  +  l~v(l)  =  0, 

t/(0)  -  l~v( 0)  =  0,  u'(l)  +  l+v(  1)  =  0, 

t/(0)  —  l~v(0)  =  0,  u'(l)  +  l~v(l)  =  0. 

We  are  interested  in  non-vanishing  solutions,  because  solutions  of  the  stationary  Burgers’ 

(%\ 

equation  are  given  by  — 2e . v  -f,  and,  therefore,  the  existence  of  a  zero  of  v(x)  together  with 

v(x) 

uniqueness  of  solutions  of  the  initial  value  problem  would  imply  a  blow-up  solution. 

We  next  introduce  the  classical  Priifer  transformation  for  (1.2.14)-(1.2.15),  namely, 

v(x)  =  r(x)  coscf)(x),  (1.2.16) 

v\x)  =  r(x)  sincf)(x).  (1.2.17) 

Differentiating  (1.2.16)  yields 

v'(x )  =  r'(x)  cos^(a;)  —  r(x )  siruf)(x)  (f>'{x).  (1.2.18) 

Now  combine  (1.2.17)  and  (1.2.18),  to  get 

r(x)sixuf)(x)  =  r'(x)  cos (f){x)  -  r{x)  sin <fi(x)  (j)'(x).  (1.2.19) 

Prom  (1.2.14)  we  have 

v"{x )  =  ( F(x )  +  \)r(x)  coscf)(x),  (1.2.20) 

and  differentiating  (1.2.17)  gives 

v"(x)=r'(x)  siiuf>(x)  +  r(x)  cos 4>(x)  <fr'(x).  (1.2.21) 

Combining  (1.2.20)  and  (1.2.21)  gives 

( F(x )  +  A )r(x)  cos <p(x)  =  r'(x)  sin <j)(x)  +  r(x)  cos cf>(x)  (j)'{x).  (1.2.22) 

Multiplying  (1.2.19)  by  cos <j>(x)  and  adding  to  (1.2.22)  times  sin0(a;)  gives 

r'(x)  =  (1  +  F(x)  -I-  A)r(x)  sin<fi(x)  cos (f>(x).  (1.2.23) 

Equation  (1.2.23)  can  be  solved  in  terms  of  0  as 

r(x)  =  r(0)  exp  [  J  (1  +  F(s )  +  A)  sin^(s)  cos (p(s)  ds).  (1.2.24) 

Multiplying  (1.2.22)  by  cos(p(x)  and  subtracting  (1.2.19)  times  sin</>(a;)  gives 

r(x)<j)'(x)  —  (F(x)  +  X)r(x)  cos 2<t>{x)  -  r(x)  si n2(f>(x).  (1.2.25) 
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Since  we  are  looking  for  solutions  u(x)  which  have  no  zeros  on  [0, 1]  and  also,  since  (1.2.24)- 
(1.2.25)  are  linear  in  r(x),  we  may  assume,  without  loss  of  generality,  that  r(x)  >  0  for 
x  G  [0,1]. 

Using  this  assumption,  we  can  divide  by  r(x)  in  (1.2.25)  to  obtain 

<jf(x)  =  (F( x)  +  A)  cos2<p(x )  —  sin^x).  (1.2.26) 

In  what  follows,  equation  (1.2.26)  is  the  main  equation  used  in  our  analysis.  On  one  hand 
it  is  independent  of  r(x ),  on  the  other  hand  r(x)  is  completely  determined  by  r(0)  and  <p{x) 
in  (1.2.24). 

Rewriting  the  boundary  conditions  (1.2.15),  taking  into  account  that  r(0)  and  r(l)  are 
not  zero,  and  after  some  simplification  using  (1.2.26),  we  obtain  the  following  four  sets  of 
boundary  value  problems 


(j)'(x) 

=  ( F(x )  +  A)  cos2<p(x)  —  sin  2<f>(x), 

(1.2.27) 

m 

=  arctan(/Sl), 

(1.2.28) 

=  —  arctan(/S2), 

(1.2.29) 

where  Si,  S2  €  {+,  -}. 

Our  main  result  is  the  following  theorem. 

Theorem  1.2.1.  The  following  statements  hold: 

a)  for  any  k  >  0  there  exists  an  antisymmetric  solution  of  the  stationary  Burgers’  system 
(1.2.3).  This  solution  is  the  unique  solution  of  (1.2.3)  for  sufficiently  large  k; 

b)  for  sufficiently  small  k,  in  addition  to  the  antisymmetric  stationary  solution,  there  also 
exist  at  least  two  non- antisymmetric  stationary  solutions; 

c)  the  antisymmetric  stationary  solution  is  asymptotically  stable  for  sufficiently  large  k. 

1.3  Harmonic  Forcing  for  Linear  Distributed  Parameter  Systems 

In  general,  the  ability  to  systematically  control  or  influence  nonlinear  effects  would  make 
a  substantial  contribution  to  existing  and  emerging  commercial  and  defense  research  and 
development  programs.  Notable  examples,  widely  appreciated  within  the  aerospace  industry, 
include  the  development  of  flight  controllers  for  high  angle-of-attack  or  high  agility  aircraft. 
Indeed,  the  importance  of  including  the  nonlinear  behavior  of  aerodynamic  parameters,  such 
as  the  coefficient  of  lift,  as  a  function  of  the  angle-of-attack  has  long  been  recognized  since 
at  high  angles-of-attack,  wind  angle  moments  also  exhibit  nonlinear  effects  which  cannot  be 
ignored.  Another  area  of  interest  is  the  control  of  flutter,  which  can  shorten  the  life  cycle  of 
aircraft  and  aircraft  parts.  Indeed,  one  example  of  the  potential  impact  of  nonlinear  control 
in  problems  of  flow  control  is  in  the  control  of  instabilities  in  the  unsteady  separated  shear 
layer,  which  has  been  experimentally  shown  to  greatly  influence  stall  and  lift  behavior  at 
high  angles  of  attack. 
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It  is  worth  noting,  however,  that  the  active  control  in  experiments,  such  as  Batill  and 
Mueller  [100],  is  based  on  a  priori  harmonic  forcing,  while  in  nonlinear  systems  with  resonance 
it  is  known  that  simple  harmonic  forcing  will  not  necessarily  produce  the  desired  response 
(see  e.g.  [41]  for  an  analysis  of  the  steady  state  response  of  nonlinear  systems  by  center 
manifold  methods).  This  point  was  also  illustrated  theoretically  by  Keefe  [59],  who  showed 
that  the  success  of  a  priori  control  for  the  Ginzburg-Landau  equation  was  dependent  on  the 
initial  state  of  the  system,  which  of  course  may  not  be  controlled  or  even  known,  and  that 
therefore  undesirable  responses  are  to  be  expected  in  the  nonlinear  regime.  Moreover,  the 
computation  by  Fuglsang  and  Cain  [38]  of  flow  over  an  open  cavity  suggests  that  harmonic 
forcing  at  non-resonant  harmonic  frequencies  can  produce  a  limit  cycle  or  chaotic  response 
that  is  far  more  severe  than  the  natural  harmonic  resonance. 

We  consider  a  special  class  of  Single  Input  Single  Output  (SISO)  linear  distributed  pa¬ 
rameter  control  systems  in  the  form 


z  —  Az  +  bu, 

(1.3.1) 

z(  0)  =  z0, 

(1.3.2) 

y  =  cz 

(1.3.3) 

where  A  is  the  infinitesimal  generator  of  a  Co  semigroup  in  a  Hilbert  space  Z  and  b  €  £(R,  Z ), 
c  €  £(Z ,  K).  Here  C(X,  Y )  denotes  the  space  of  bounded  operators  from  X  to  Y. 

We  assume  that  the  input  u  is  given,  in  feedback  form,  as  the  output  of  a  harmonic 
oscillator  with  frequency  a: 


w  =  Sw, 


tu(0)  = 


u  -  Tw, 


(1.3.4) 


where  T  is  a  given  1x2  matrix:  T  =  [71, 72].  Thus  u  represents  a  periodic  function  of  period 
T  =  27t/q:  as  a  linear  combination  of  sin(at)  and  cos  (at),  namely, 

u(t)  =  71  sin(o!t)  +  72  cos  (at).  (1.3.5) 

Problem  1.3.1.  Suppose  that  we  are  given  the  input  u  in  (1.3.5),  find  an  initial  condition  z0 
in  (1.3.2)  so  that  the  output  y  in  (1.3.3)  is  a  nontrivial  periodic  function  of  period  T  =  27r/a. 

It  is  well  known  in  finite  dimensional  linear  control  theory  that  if  a  system  is  driven  by  a 
periodic  input  for  which  the  complex  frequency  ia  is  a  transmission  zero  of  the  system,  then 
the  output  of  the  system  is  zero  for  all  time.  Therefore,  we  should  also  state  the  following 
more  general  problem. 

Problem  1.3.2.  Find  conditions  on  ( A,b,c )  guaranteeing  there  is  a  nontrivial  periodic  out¬ 
put  with  the  desired  period  for  all  a  and  arbitrary  71,  72  with  7i  +  7I  ^  0,  i.e.,  that  the 
system  will  support  a  periodic  output  of  arbitrary  period. 
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In  order  that  the  solution  to  (1.3.1)  be  periodic  we  would  at  least  need  that  z(T )  =  z(0). 
By  the  variation  of  parameters  formula  we  have 

rt 


z(t)  =  eAtz0  +  [  e A(~l  T^u(r)  dr. 
Jo 

In  order  that  z{-)  satisfy  z{T )  =  z(0)  =  z0  we  need 

(l-eAT)z0=  [  eA(T~T)bu{T)dT. 
Jo 


(1.3.6) 


(1.3.7) 


Assumption  1.3.1.  In  this  paper  we  will  avoid  the  various  technical  difficulties  and  make 
the  assumption  that  A  is  a  discrete  Riesz  spectral  operator  with  simple  eigenvalues  (  mul¬ 
tiplicity  one)  and  eigenvectors  {' ipj These  eigenvectors  form  a  Riesz  basis  in 

Z  (i.e.,  a  linear  isomorphic  image  of  an  orthonormal  basis).  In  this  case  the  adjoint  A*  is 
also  a  discrete  Riesz  spectral  operator  whose  eigenvectors  form  a  biorthogonal  Riesz 

basis,  i.e.,  (ipj,  tpl)  =  6jk. 

Note  that  due  to  our  assumption  that  b  and  c  are  bounded  rank  one  operators,  we  have 
a  well  defined  transfer  function  given  by  g(s)  =  c(sl  -  A)~1b. 


Assumption  1.3.2.  A  natural  assumption  on  our  system  is  that  the  transfer  function  is 
real,  i.e., 


g(s)  =  g(s). 


(1.3.8) 


For  systems  governed  by  differential  equations  with  real  coefficients  this  condition  is  auto¬ 
matic. 


Definition  1.3.1.  A  complex  number  sq  is  a  transmission  zero  if  g(s0)  =  0. 

Assumption  1.3.3.  Our  final  assumption  is  that  there  are  no  pole  zero  cancellations.  That 
is,  we  assume  that  if  s0  is  a  transmission  zero,  then  so  6  p{A),  the  resolvent  set  of  A. 

Our  main  results  are  stated  in  the  following  theorem. 

Theorem  1.3.1.  Let  the  operator  A  in  (1.3.1)  be  a  discrete  Riesz  spectral  operator  with 
a(A)  =  {Aj}^,  the  input  u  is  given  by  (1.3.5)  with  7?  +  ^  0  and  let  ( A,b,c )  satisfy 

Assumptions  1.1,  1.2  and  1.3.  Then  we  have  the  following  results. 

1.  There  exists  an  initial  condition  zq  so  that  the  solution  z  to  (1.3.1)  is  periodic  with 
period  T  —  2r:lo.  provided 

distance(a(A),  {kai  |  k  —  0, ±1, ±2,  •••})>0.  (1.3.9) 


Furthermore,  the  system  supports  all  positive  periods  T  (i.e.,  we  can  find  a  periodic 
solution  for  all  possible  frequencies  a)  if 

distance{a{A),^f)  >  0 

where  C°  =  {A  €  C  :  Re  A  =  0}  denotes  the  imaginary  axis. 
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2.  In  this  case,  there  is  a  nontrivial  periodic  output  y  if  and  only  ifia  is  not  a  transmission 
zero,  i.e.,  g(ia)  7^  0. 

3.  Finally,  let  us  denote  the  amplitude  of  the  periodic  input  u  by 

Mu  =  sup  |u(t)|  =  + 

t€[0,T]  v 

Then  the  amplitude  of  the  output  y  is  a  linear  function  of  the  amplitude  of  the  input 
u.  In  particular,  the  output  can  be  written  in  the  forms 

y(t)  =  [Res(ia)]u(£)  +  ^[lmg(ia)]^(t)  (1.3.10) 

=  Mu\g(ia)\[ii  sin(af)  +  72 cos(a£)] 

=  Mu\g(ia)\ sin(of£  +  <j>)  (1.3.11) 

where  71 2  +  72  2  =  1  and  we  can  easily  write  explicit  formulas  for  fi,  72  and  (p  in  terms 
of  71,  72  and  g(ia).  Thus  the  amplitude  My  of  y  can  be  written  as 

My=  sup  \y(t)\  =  Mu\g(ia)\. 
te[o,r] 

The  proof  is  based  on  a  functional  calculus  valid  for  the  discrete  Riesz  spectral  operators 
considered  here.  In  particular,  under  the  assumptions  of  the  theorem,  we  can  use  spectral 
theory  to  obtain  explicit  representations  for  the  initial  data  z0  and  the  solution  z.  We  show 
that 


z0  =  ~('7ia  +  j2A){A2 +  a2)  lb, 


(1.3.12) 


z{t)  = 


sin(o:£)  (—71 A  +  72a)  +  cos(a£)  (—71a  -  72^) 


(. A2  +  a2)~1b . 


(1.3.13) 


Applying  c  to  (1.3.13),  using  the  resolvent  identity  and  the  fact  that  the  transfer  function 
is  real  (cf,  (1.3.8)),  we  can  write  y  as 

y(t)  =  sin(crt)  [71  Re^(ia)  -  72  Im^(*a)] 

+  cos  (at)  [71  Im  g{ia)  +  72  Rep(za)]  (1.3.14) 

from  which  a  straightforward  calculation  shows  that  y  is  nontrival  if  and  only  if 

\g(ia)\2  =  (R eg(ia))2  +  (lmg{ia))2  ±  0, 

i.e.,  if  and  only  if  ia  is  not  a  transmission  zero. 

Part  3  of  the  theorem  follows  from  formula  (1.3.14).  In  particular,  a  direct  calculation 
shows  that  the  applitude  My  of  y  is  related  to  the  amplitude  Mu  of  u  by 

My  =  (7i  +  7 2)  (Re^'a)2  +  Im^fa)2) 

=  (7i  +7l)  bMI2  =  Ml\g{ia) |2. 
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1.4  Existence  of  Numerical  Equilibria 

Another  result  obtained  during  this  research  effort  was  concerned  with  the  effect  of  boundary 
control  on  the  structure  of  attractors  for  Burgers’  equation.  We  were  able  to  provide  a  rig¬ 
orous  proof  that  the  structure  of  the  attractor  is,  in  general,  nontrivial.  Indeed,  for  a  special 
class  of  forcing  terms  we  are  able  to  show  that  as  certain  gain  parameters  in  the  boundary 
feedback  control  are  increased  from  small  positive  values  to  large  positive  values,  the  number 
of  stationary  solutions  vary  from  three  to  one.  This  result  supports  our  conjecture  that  for 
large  values  of  the  gains,  there  is  a  single  global  asymptotically  stable  equilibrium.  In  par¬ 
ticular,  for  zero  external  forcing  term  this  would  say  that  the  resulting  closed  loop  system 
is  asymptotically  stable.  This  conjecture  is  based  on  development  of  the  notion  of  nonlinear 
zeros  as  a  nonlinear  enhancement  of  the  classical  concept  of  rootlocus.  In  particular,  we  have 
given  a  definition  of  the  zero  dynamics  associated  with  the  closed  loop  Bugers’  system  and 
shown  that  this  zero  dynamics  systems  has  a  single,  global,  asymptotically  stable,  equilib¬ 
rium.  One  would  therefore  expect  that  as  the  gains  are  increased,  trajectories  of  the  closed 
loop  system  would  approach  the  corresponding  trajectories  of  the  zero  dynamics. 

During  this  period,  in  joint  work  with  Dr.  John  Burns  (VT  &  State  University),  we  also 
made  an  interesting  discovery  concerning  the  long  time  behavior  of  solutions  to  Burgers’ 
equation  on  a  finite  interval  with  Neumann  boundary  conditions.  It  is  easy  to  see  that,  for 
this  problem,  constants  are  equilibria,  and  for  the  related  linearization  about  zero  -  the  one 
dimensional  heat  equation  with  Neumann  boundary  conditions  -  it  is  well  known  that  the 
steady  state  temperature  is  a  constant.  Namely,  the  steady  state  temperature  distribution  is 
the  mean  value  of  the  initial  temperature  distribution.  For  the  Burgers’  equation  with  small 
initial  data  this  same  type  of  result  holds  as  a  consequence  of  the  Center  Manifold  Theorem. 
That  is,  Burgers’  equation  with  Neumann  boundary  conditions  possesses  a  one  dimensional 
center  manifold  (constants)  and  it  can  be  shown  that  for  a  small  initial  conditions  solutions 
converges  exponentially  to  constant  values.  As  we  have  shown  in  our  earlier  work,  in  contrast 
to  the  heat  equation,  the  steady  state  constant  is  not  simply  the  mean  of  the  initial  condition, 
but  it  also  depends  in  some  complicated  way  on  both  the  viscosity  parameter  and  the  shape 
of  the  initial  condition. 

Since  the  Center  Manifold  Theorem  is  only  a  local  result,  a  natural  question  is  whether, 
for  arbitrary  initial  data,  the  corresponding  solution  of  Burgers’  equation  tends  to  a  constant 
steady  state.  The  answer  to  this  question  is  still  unresolved.  Nevertheless,  we  were  able  to 
resolve  in  the  affirmative  an  intermediate  question.  If,  for  a  given  initial  condition,  the 
solution  approaches  a  time  independent  steady  state,  even  in  the  L  sense,  then  this  steady 
state  must  be  a  constant. 

In  spite  of  this  result,  after  considerable  numerical  testing,  it  was  discovered  by  researchers 
at  Virginia  Tech  University,  that  for  moderately  small  viscosity  and  larger  “antisymmetric” 
initial  conditions,  numerical  solutions  can  approach  a  nonconstant,  time  independent  steady 
state.  We  could  only  conclude  that  this  is  simply  a  numerical  anomaly.  Nevertheless,  due 
to  the  relevance  of  hydrodynamic  problems  in  applications,  it  was  important  to  understand 
how  these  nonconstant  numerical  stationary  solutions  arise.  Our  discovery  was  that  these 
numerical  solutions,  for  a  fixed  mesh  size  or  degree  of  approximation,  approach  explicit 
solutions  of  the  equation  that  satisfy  the  boundary  conditions  only  to  within  values  that 
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are  approximately  machine  precision  zero  (or  smaller).  In  particular,  these  numerical  sta¬ 
tionary  solutions  are  given  in  terms  of  hyperbolic  tangent  functions  whose  derivatives  at 
the  end  points  behave  like  hyperbolic  secant  squared.  For  moderate  initial  conditions  these 
derivatives  are  on  the  order  of  10-30  which  are  well  below  machine  precision  zero. 

This  work  sheds  important  light  on  the  use  of  numerical  proofs.  In  particular,  numerical 
calculations  have  been  used  to  “suggest”  that  Euler  equations  do  not  have  unique  solutions. 
The  justification  for  this  claim  is  that  a  “very  fine  mesh”  is  used  in  the  calculation.  The 
conclusion  of  our  work  is  that  “numerical  based”  proofs  of  non-existence  must  be  done  with 
extreme  care. 

Burgers’  equation  on  the  interval  (0, 1)  subject  to  Neumann  Boundary  Condition  is  given 
by  the  dynamical  system 

wt  -  ewxx  +  wwx  =  0,  (1.4.1) 

x  £  (0, 1),  t  >  0 
wx(0,t)  =  wx(l,t)  =  0 
w(x,  0)  =  (j)(x ), 

The  associated  stationary  Burgers’  problem  is 

-evxx  +  vvx  =  0,  (1.4.2) 

v*(0)  =  vx(l)  =  0. 

It  is  not  clear  without  further  information  that  solutions  of  (1.4.1)  should  even  exist  for  all 
time.  The  answer  to  this  question  is  contained  in  the  recent  work  [119]  which,  for  the  special 
case  of  (1.4.1),  gives  the  following  result. 

Theorem  1.4.1.  [119]  For  arbitrary  initial  data  <p  £  L2( 0, 1)  and  0  <  T  <  oo, 

a)  (1-4-1)  has  a  unique  weak  solution 

w  £  L°°([ 0,  T],  L2(0, 1))  n  L2([0,  T],  h\ 0, 1)), 

b)  On  any  cylinder  [0, 1]  x  [t0,T]  for  any  0  <  to  <  T  <  00,  w  £  if2,1([0, 1]  x  [t0,T]). 

c)  There  is  a  globally  defined  dynamical  system  on  the  state  space  L2{ 0, 1)  given  in  terms 
of  a  nonlinear  semigroup  {Tt,t  >  0}. 

*  Tt  is  continuous  in  t  and  £  L2{ 0, 1). 

*  Tt  is  compact  for  t  >  0. 

*  The  system  is  globally  Lyapunov  stable. 

*  There  is  a  global,  locally  compact  attractor. 

Even  though  we  cannot,  at  this  time,  say  that  solutions  of  (1.4.1)  approach  a  constant 
steady  state,  we  prove  the  following  intermediate  result. 


15 


Theorem  1.4.2.  Fix  e  >  0  and  (f>  €  L2( 0, 1).  Let  w(-,t)  be  a  weak  solution  of  (1.4-1)-  If 
there  is  a  function  h  €  L2(0, 1)  such  that 

lim  ||iu(-,£)  —  /i(-)||l2(o,i)  ~*t—> 


Then  h(-)  =  c ^  for  some  constant  c^. 

For  fixed  e  and  for  small  initial  data  numerical  approximation  of  the  solutions  to  (1.4.1) 
supports  the  conclusion  of  the  Center  Manifold  Theorem,  namely,  solutions  tend  to  a  constant 
as  t  tends  to  infinity.  But  for  small  e  and  “certain”  initial  data  (not  to  small),  the  numerical 
solution  converges  to  a  nonconstant  function,  cf.  [94].  We  are  lead  to  conjecture  the  existence 
of  some  type  of  Numerical  Stationary  Solutions  for  the  problem  (1.4.1). 

One  class  of  initial  data  for  which  we  obtain  this  anomaly  are  functions  in  the  class  S 
consisting  of  “antisymmetric”  functions,  that  is, 

S={<f)E  L2{ 0, 1)  :  <p{x)  =  -4>{  1  -  x)}  .  (1.4.3) 


For  initial  data  (f>  E  S,  a  straightforward  consequence  of  Theorem  2.1,  is  that  w(-,t)  E  S 
for  all  t  and  hence  w(l/2,t)  =  0  for  all  t  >  0.  Thus  if  limu;(a;,  t)  =  c^)£  exists  then  the 

constant  c ^  must  be  zero. 

The  nonconstant  solutions  in  S  of  the  stationary  Burgers’  equation  (not  the  boundary 
conditions)  are  given  in  terms  of  a  one  parameter  family  depending  on  the  parameter  Cq. 


h(x)  =  y/2co  tanh  (1/2  -  a:)^  ,  h  £  S. 


(1.4.4) 


It  is  easy  to  see  that  for  suitable  initial  data  and  Co  the  functions  (1.4.4)  are  actually 
numerical  stationary  solutions  to  (1.4.1),  i.e.,  they  satisfy  the  Burgers’  equation  and  they 
approximately  satisfy  the  boundary  conditions. 

Namely,  the  functions  in  (1.4.4)  satisfy  (1.4.2)  and  the  boundary  conditions 


h'{  0)  =  h'(  1)  =  — —  sech2 


=  “7, 


(1.4.5) 


where  7  is  an  exponentially  small  positive  number. 

There  is  no  reason  to  believe  that  numerical  solutions  to  Burgers’  equation  should  ap¬ 
proach  a  function  of  the  type  (1.4.4),  especially  in  light  of  Theorem  2.2  which  suggests  they 
should  approach  a  constant.  Nevertheless,  this  does  happen  for  larger  initial  data  and/or 
smaller  e. 

For  a  small  positive  number  7  there  are  exactly  two  solutions  (1.4.4)  satisfying  (1.4.5), 
i.e.,  there  exist  cj  «  0  and  (%  »  0  giving 


hj(x)  = 


2cq  tanh 


(1.4.6) 
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=  7- 


(1.4.7) 


The  solution  hx  is  very  nearly  the  zero  function,  whereas  the  solution  h?  is  not  usually  small. 
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S.  Gilliam 


October  1995: 

“Robust  Regulation  of  Nonlinear  Systems,”  Invited  talk  presented  by  Prof.  C.I.  Byrnes 
at  VPI  &  State  University. 

“Robust  Regulation  of  Nonlinear  Systems,”  Plenary  talk  presented  by  Prof.  C.I.  Byrnes 
at  NCSU  Southeast  Regional  Differential  Equations  Conference. 

“On  the  Dynamics  of  Boundary  Controlled  Convective  Reaction  Diffusion  Equations,” 
Invited  talk  presented  by  D.S.  Gilliam  in  the  minisymposium  “Control  and  Compu- 
tatuional  Fluid  Dynamics”  at  the  1995  SIAM  Annual  Meeting  at  Charlotte,  North 
Carolina. 

“Electromagnetic  Waves  Over  a  Three  Dimensional  Conducting  Half  Space  Invited 
lecture  presented  by  D.S.  Gilliam  in  the  Center  For  Research  in  Scientific  Computa¬ 
tion  at  North  Carolina  State  University. 


November  1995: 

“Robust  Regulation  of  Nonlinear  Systems,”  Invited  talk  presented  by  Prof.  C.I.  Byrnes 
in  the  Departemnt  of  Electrical  Engineering  at  the  University  of  Hong  Kong,  Hong 
Kong. 
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“Output  Regulation  Revisited”  Distinguished  lectures  series,  Department  of  Electrical 
Engineering,  IOWA  State  University,  presented  by  Dr.  Alberto  Isidori. 


December  1995: 

“Disturbance  Attenuation  for  a  Class  of  Nonlinear  Non-Minumum  Phase  Systems,” 
Invited  lecture  presented  by  Prof.  Alberto  Isidori  at  the  34th  Control  and  Decision 
Conference. 

“High  Gain  Limit  for  Boundary  Controlled  Convective  Reaction  Diffusion  Equations” 
lecture  presented  by  D.S.  Gilliam  at  the  34th  Control  and  Decision  Conference. 

“Well-Posedness  for  Controlled  Nonlinear  Damped  Membranes  with  Fixed  Boundary” 
lecture  presented  by  D.S.  Gilliam  at  the  34th  Control  and  Decision  Conference. 

January  1996: 

“The  Regulator  Problem  for  Linear  Distributed  Parameter  Systems,”  Invited  lecture 
presented  by  D.S.  Gilliam  in  the  Department  of  Systems  Science  and  Math  at  Wash¬ 
ington  University  in  St.  Louis. 

March  1996: 

“Complete  Parameterization  of  Solutions  to  the  Rational  Covariance  Problem,”  Invited 
lecture  presented  by  C.I.  Byrnes  at  Ecole  Normale  Superieure,  Paris. 

April  1996: 

“Complete  Parameterization  of  Solutions  to  the  Rational  Covariance  Problem,”  Ple¬ 
nary  lecture  presented  by  C.I.  Byrnes  at  Conference  on  the  Future  Directions  in  Applied 
Mathematics  at  Notre  Dame. 

June  1996: 

“The  Structure  of  Attractors  for  a  Boundary  Controlled  Viscous  Burgers’  Equation,” 
Invited  talk  presented  by  Andras  Balogh  at  the  96  Mathematical  Theory  of  Networks 
and  Systems  (MTNS)  Conference,  St.  Louis,  MO. 

“Recent  Results  on  Boundary  Control  for  a  Class  of  Higher  Dimensional  Convection- 
Reaction  Diffusion  Equations,”  Invited  talk  presented  by  David  S.  Gilliam  at  the  96 
Mathematical  Theory  of  Networks  and  Systems  (MTNS)  Conference,  St.  Louis,  MO. 

“Robust  Nonlinear  Output  Regulation,”  Invited  talk  presented  by  Prof.  Alberto  Isidori 
at  the  96  Mathematical  Theory  of  Networks  and  Systems  (MTNS)  Conference,  St. 
Louis,  MO. 

August  1996: 

“On  the  Rational  Caratheodory  Extension  Problem,”  Invited  lextrue  presented  by 
Prof.  C.I.  Byrnes  at  the  Fifth  Bozeman  Conference  on  Computation  and  Control. 
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“On  the  Dynamics  and  Control  of  a  Class  of  Higher  Dimensional  Convection  Reac¬ 
tion  Diffusion  Equations,”  Invited  lecture  presented  by  David  Gilliam  at  the  Fifth 
Bozeman  Conference  on  Computation  and  Control. 

“Numerical  Stationary  Solutions  for  Burgers’  Equation  with  Neumann  Boundary  Con¬ 
ditions,”  Invited  lecture  presented  by  David  Gilliam  at  the  Fifth  Bozeman  Conference 
on  Computation  and  Control. 

October  1996: 

“Zeros  and  Zero  Dynamics  for  Linear  Distributed  Parameter  Systems,”  A  series  of  lec¬ 
tures  given  by  Dr.  David  Gilliam  in  the  Texas  Tech  University  Joint  Civil  Engineering, 
Mathematics  and  Mechanical  Engineering  Seminar. 

December  1996: 

“Global  normal  forms  for  MIMO  nonlinear  systems,  with  application  to  stabilization 
and  disturbance  attenuation”  Lecture  presented  by  Dr.  Alberto  Isidori  at  the  35rd 
IEEE  Conf.  Decision  and  Control  (Kobe,  Japan). 

“Z/2  disturbance  attenuation  and  performance  bounds  for  linear  non-minimum  phase 
square  invertible  systems”  Lecture  presented  by  Dr.  Alberto  Isidori  at  the  35rd  IEEE 
Conf.  Decision  and  Control  (Kobe,  Japan). 

“Immersions  and  the  internal  model  principle:  tools  for  robust  nonlinear  control”  Lec¬ 
ture  presented  by  Dr.  Christopher  I.  Byrnes  at  the  35rd  IEEE  Conf.  Decision  and 
Control  (Kobe,  Japan). 

“The  rational  covariance  extension  problem  with  applications  to  speech  synthesis” 
Invited  lecture  presented  by  Dr.  Christopher  I.  Byrnes  at  The  University  of  Tokyo. 

“Robust  nonlinear  control”  Invited  lecture  presented  by  Dr.  Christopher  I.  Byrnes 
at  Dewan  Riset  Nasional,  Serpong,  Indonesia. 

January  1997: 

“Global  L2~ Gain  State  Feedback  Design  for  a  Class  of  Nonlinear  Systems”  Lecture 
presented  by  Dr.  Alberto  Isidori  at  COSY  Workshop  on  Nonlinear  Systems  held  at 
the  ETH-Zurich,  Switzerland. 

February  1997: 

“The  geometry  of  positive  real  functions  with  applications  to  signal  processing  and  to 
speech  synthesis”  Invited  lecture  presented  by  Dr.  Christopher  I.  Byrnes  at  The  Texas 
Tech  University. 

March  1997: 

“The  rational  covariance  extension  problem  with  applications  to  speech  synthesis”  In¬ 
vited  lecture  presented  by  Dr.  Christopher  I.  Byrnes  at  Texas  Tech  University. 
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“Robust  Output  Regulation  of  Nonlinear  Systems”  Invited  lecture  presented  by  Dr. 
Alberto  Isidori  at  the  University  of  Aalborg,  Denmark. 

April  1997: 

“On  the  Global  Dynamics  of  a  Class  Boundary  Controlled  Hydrodynamic  Systems  in 
Higher  Dimensions,”  Lecture  presented  by  Dr.  Victor  I.  Shubov  at  the  Texas  Partial 
Differential  Equations  Conference  in  Denton,  TX. 

May  1997: 

“Nonlinear  control  systems”  Invited  lecture  presented  by  Dr.  Christopher  I.  Byrnes  at 
the  1997  AFOSR  Contractors  Review  Meeting,  Wright  Patterson  AFB,  Dayton,  OH. 

June  1997: 

“On  the  covariance  extension  problem”  5th  IEEE  Mediterranian  Conference,  Paphos, 
Cyprus  Invited  lecture  presented  by  Professor  C.  I.  Byrnes 

“Robust  Output  Regulation  of  Nonlinear  Systems”  Invited  lecture  presented  by  Dr. 
Alberto  Isidori  at  the  Imperial  College,  London,  UK. 

“Robust  Nonlinear  Control”  Invited  lecture  presented  by  Dr.  Christopher  I.  Byrnes  at 
the  28th  AIAA  Fluid  Dynamics  Conference/  4th  AIAA  Shear  Flow  Control  Conference, 
Snowmass  Village,  CO 

July  1997: 

Dr.  David  S.  Gilliam  participated  in  12th  International  Conference  of  the  American 
Meterological  Society  on  Boundary  Layers  and  Turbulence,  Vancouver,  BC. 

“On  Duality  between  Filtering  and  Interpolation”  University  of  Tokyo,  Japan  Invited 
lecture  presented  by  Professor  Anders  Lindquist. 
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